Quantum heat engines and nonequilibrium temperature 



Ramandeep S . Johal 
Department of Physics, Indian Institute of Science Education and Research Mohali, 
Transit Campus: MGSIPA Complex, Sector 26, Chandigarh 160019, India 
(Dated: September 24, 2009) 

A pair of two-level systems initially prepared in different thermal states and coupled to an external 
reversible work source, do not in general reach a common temperature at the end of a unitary work 
extraction process. We define an effective temperature for the final nonequilibrium but passive state 
of the bipartite quantum system and analyse its properties. 
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I. INTRODUCTION 

Consider the thermodynamic problem of work extrac- 
tion [l|, |2| from two systems at different temperatures 
Ti and T2 (let Ti > T2) by coupling them with a re- 
versible work source. It is assumed that internal energy 
of each system is Ui = CiTi, where Ci is independent of 
temperature. The process of work extraction stops when 
the two systems reach a common final temperature Tf. 
Work performed is given by the difference of initial and 
final energies: 



Wo = CiTi + C2T2 - {Ci + C2)Tf. 



(1) 



Now to extract maximal work, the process is assumed 
to be thermally isolated in which thermodynamic en- 
tropy of the total bipartite system is preserved. This 
criterion yields the value of the final temperature as 
Tf = (ri)«/(i+«)(T2)i/(i+«), where ^ = C1/C2.. 

One can discuss a cyclic process which proceeds in the 
following two steps: i) the two systems prepared as above 
and coupled to a reversible work source, are used to ex- 
tract an amount of work given in Eq. ([1]) whence the 
systems reach a common final temperature; ii) the sys- 
tems are then brought back to their initial states by sep- 
arating them from the work source and making contact 
with thermal baths at Ti and T2 respectively. In the 
second step, the system 1 absorbs heat from the hotter 
bath and system 2 rejects some heat to the cold bath. 
The efficiency of this cyclic process is 



16- 



(2) 



where 6 = T2/T1. This system behaves very similar to a 
cycle discussed by Leff 3] which is made up of a sequence 
of infinitesimal Carnot cycles and where both the heat 
source and the sink have a finite heat capacity. 

The problem of work extraction has also been ad- 
dressed from a quantum mechanical point of view 0, 
0]. Although the possibility of a quantum heat engine 
and validity of thermodynamic bounds has been recog- 
nised since 1950s Q, the recent developments in nan- 
otechnology and quantum information processing have 
contributed to enhanced interest in quantum thermody- 
namic machines [1, [l^l- Alongside, such models provide 



insight into fundamental questions about thermodynam- 
ics such as Maxwell's demon and universality of the sec- 
ond law [1, [111, [I2| • Many models employ few- level quan- 
tum systems as the working medium, such as quantum 
harmonic oscillators, spin-systems, particle-in-box and so 
on [H, [3, [1^ [H, [l3| . Usually the cycle is a quantum 
generalization of the well-known classical Carnot, Otto, 
Brayton heat engines which follow four-step cycles. In 
another class of models, instead of the two classical or 
macroscopic systems as discussed in the preceding para- 
graph, one can form a two-step engine using two quantum 
systems (6|. Recently, such a quantum heat engine em- 
ploying two two-level systems (TLS) was discussed and 
implications of the optimization of work on the structure 
of the engine were highlighted [Isj . 

However, quantum engines being small systems, the 
validity of thermodynamic behaviour is not guaranteed. 
For instance, after work extraction in the latter class 
of models, the two systems may not reach mutual equi- 
librium. In this paper, we further discuss the two-step 
model for work extraction using two TLS, focusing on the 
final passive state (which is in general a nonequilibrium 
state) from a thermodynamic perspective. We define an 
effective temperature for this state and analyse its prop- 
erties. The paper is organised as follows. In section II, 
we introduce the model of quantum heat engine. In sec- 
tion IIA, the temperatures of subsystems are evaluated; 
the validity of thermodynamic definitions is enforced by 
deriving the specific heats of subsystems in section IIB. 
Section III proposes a definition for effective temperature 
of the composite system, which is calculated explicitly in 
different regimes of parameter values. We also compare 
some of the other definitions in literature for effective 
temperature of nonequilibrium systems, in section IV. 
Concluding ideas are given in section V. 



II. QUANTUM MODEL FOR WORK 
EXTRACTION 



Consider two TLS labeled R and S with hamiltonians 
Hn and Hs, prepared in thermal states pn and ps cor- 
responding to temperatures Ti and T2. The hamiltonian 
of the total system is H = ^ I + I (E) Hs- The ini- 
tial state of the composite system is pin = Pr® Ps- The 
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eigenvalues of H are {0, 02, ai, oi +02} given that energy 
eigenvalues of Hji and Hs are (0, Oi) and (0, 02), respec- 
tively. The eigenvalues of the initial density matrix are 
{risi, riS2, f2Si, ^282}. Here the probability to find each 
system in its excited state is 



1 



r2 



1 + e'^i/'^i ' 



S2 



1 



/T2 



(3) 



with ground state probabilities being ri = (1 — r2) and 
si — {1 — 82)- We set Boltzmann constant fee — 1. 

The initial mean energy of the composite system is 
U = air2 + a2S2- Let us for concreteness choose, ai > a2. 
Within the approach based on quantum thermodynam- 
ics, the process of work extraction is a unitary process 
which preserves not only the magnitude of the entropy 
but also all eigenvalues of the density matrix describ- 
ing the state of the system. It has been shown in earlier 
works that under such a process, the state which 

corresponds to a minimum value of the final energy is 
Pj = PS ® PR, with eigenvalues {riSi, r2Si, riS2, r2S2}. 
Effectively, it means that the two systems exchange or 
swap their initial probability distributions in the final 
state. In other words, work performed is maximum if 
U = aiS2 + a2r2 and is given by 



>V(ai,a2) = C/ - U ={ai-a2){s2-r2). 



(4) 



Net work is extracted if W < which requires the fol- 
lowing condition: 



S2 < r2 



T2 02 



(5) 



The efficiency of this engine is 77 = I — j which is inde- 
pendent of temperature and its upper bound is Carnot 
value. 



A. Temperatures of subsystems after work 
extraction 

Now we study temperatures in the final state. After 
work, the mean energy of subsystem l\sUi= aiS2, and 
U2 — (12^2 ■ Let us consider two such set-ups specified by 
the pair of energy parameters (ai , 02) and (oi -I- dai , 02 -I- 
da2). Comparing the final states after work extraction, 
the change in energy of subsystem I is 



' cls2 

dUi — S2dai + ai- — da2- 
da2 



(6) 



We follow the standard interpretation of work as the 
change in mean energy due to shift in energy levels, at 
constant probabilities [H, HH . Similarly, heat is de- 
fined to be the change in mean energy when the energy 
levels stay fixed, but probability of occupation changes. 
Thus the heat contribution for system I is given by 



dQ, 



ds2 , 
ai- — da2- 
da2 



(7) 



Similarly for subsystem 2, we have 



dQ2 



dr2 
0.2-, — dai. 
da I 



(8) 



Let us now study entropy of each subsystem. In the 
initial state, the entropy of subsystems are given by 
Si = — (ri Inri + r2 lnr2) and 6*2 = — (si Insi + S2 lns2) 
respectively. After work, due to exchange of probabili- 
ties between the subsystems, we have = S2, S2 — Si. 
Thus for subsystem, say I, the change in entropy of the 
final state under a variation of the parameter 02 is 



dSi = dS2 

a2 ds2 



da2- 



T2 da2 

Now we evaluate the final temperature of system 1 as 



(9) 
(10) 



dQi rj. ai 

r — ^ 2 — • 

dS^ a2 



Similarly, we obtain for system 2 



±2 — il . 

a-i 



(11) 



(12) 



These values of temperatures are precisely which may be 
obtained directly from the final probability distributions 
of the TLS, because a TLS can always be assigned an 
effective temperature. 

Using Eq. ([5|), it can be seen that after work ex- 
traction, the hotter subsystem I cools down {T^ < Ti), 
where as the relatively cold subsystem 2 now has a higher 
temperature (Tj > T2). Note that the sign of differ- 
ence {Ti —T2) is not determined; it is possible to have 
(Tj^ < T2). But this does not violate the second law, 
because the condition (O also ensures that energy flows 
from the hot to the cold system. Thus the change in 
energy of system 1, AUi = ai(s2 — 7*2) < and the cor- 
responding change in system 2 is /S.U2 = ai(r2 — S2) > 0. 



B. Heat capacity of subsystems 

The canonical heat capacity of subsystem 1 in the fi- 
nal state is related to the fluctuations of energy in a 
well-known way S]. However, heat capacity may also 
be evaluated as follows. Consider the final temperature 
as function of oi and a2 (Eq. (jlip ). Then a change in 
temperature resulting from a variation in these parame- 
ters is 



T2 aiT2 
dT^ = — dai ~ 7 — 
0.2 (02) 



da-2 



(13) 



Then keeping ai fixed (which is equivalent to keeping 
volume of subsystem 1 fixed, because change in oi for 
subsystem f in the final state is interpreted as work, see 
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Eq. ([6])), the heat capacity (at constant volume) in the 
final state of system 1 is 




where we have used the following identity 

C2T2 



dS2 

da2 



(14) 
(15) 

(16) 



Here C2 is the canonical heat capacity of the subsystem 
2 in its initial state at temperature T2. Similarly, we get 
the result C2 — Ci. Thus upon swap-transformation, the 
specific heats of the two subsystems also get exchanged. 

To recapitulate, the standard thermodynamic process 
in which two macroscopic bodies at different tempera- 
tures are coupled to a work source, the final tempera- 
tures of the two bodies are said to become equal. In the 
quantum framework, the subsystems in general do not 
reach mutual thermal equilibrium. In the next section, 
we ask: can the whole bipartite system be characterised 
by a global effective temperature in the final state, even 
though it is a nonequilibrium state with subsystems at 
different temperatures ? 



III. 'TEMPERATURE' FOR THE BIPARTITE 
SYSTEM 

For subsystem i, we observed in the previous section 
that temperature can be defined thermodynamically. In 
this section, we extend the thermodynamic definition to 
the nonequilibrium final state of the composite system. 

For convenience, we define 02/01 — v. So the final 
temperatures are rewritten as: = T2/1' and = Tiv. 
Thus for given reservoir temperatures (Ti ,T2)^ the final 
temperatures of subsystems is determined by a single pa- 
rameter V, which is also related to the efficiency of the 
engine v = 1 — rj. Consider different final and initial 
states which are characterized by the same parameter v, 
but which may yield different amounts of work. 

At a given value of v, the changes in ai and 02 are 
related as 



da2 — vdai. 



(17) 



Thus the heat exchanged by system 1 in such a process 
can be rewritten from Eq. ([7]) as 



dQi = ai—^vdai 
da2 

dS2 , 

— a2- — dai. 
da2 



(18) 
(19) 



Also, Eq. ([H]) yields dQ2 Then the total heat exhanged 
by bipartite system is dQ = dQi + dQ2. 



Similarly, the von-Neumann entropy of the bipartite 
system is the sum of subsystem entropies, S =81 + 82 = 
81 + 82, and the variation in total entropy is 



d8 



d8i 
da I 



dS2 
da2 



dai, 



(20) 



at a given using Eq. (fT7)) . Then we define the effective 
temperature as the ratio of heat variation to the entropy 

variation, T ee (^^^ , yieldir 



Img 



T : 



02 



\ da2 dai J 



(02)" dS2 
aiT^ da2 



ai dr2 
Ti dai 



(21) 



Using Eq. (HH) and - ^ = we finally get 



T = 



C2T[ + C1T2 

{Ci + C2) 



(22) 



The above formula is the main result of the present paper. 
It resembles the thermodynamic expression if the two 
systems at temperatures and T2 with constant heat 
capacities C2 and Ci respectively, come to a common 
temperature T, without doing any work, (see Eq. ((T])). 
Now we evaluate the effective temperature for the case 
of two TLS, and discuss its features. The canonical heat 
capacity of a TLS is given by the well known expression 



exp[Q,/Ti] 
{l + eMa^/T^W 



(23) 



We first discuss the limit when ai/Ti = x ^ 1. Then 



74 and the ratio C1/C2 = C ^ [Glv f - Thus 



the efficiency is given by 77 = 1 — The effective tem- 
perature in this regime is 



T 



(1+0 



ii + e). 



(24) 



At Carnot hmit, ^ ^ 1 and T/Ti = (1 + e)/2. How- 
ever, note that this formula holds in general also, because 
Carnot limit implies ai/Ti — > 02 /T2 and so C1/C2 — + 1. 
Here the extracted work is vanishingly small and the final 
temperature is expected to be (Ti -I- T2)/2. 

At the other extreme, for ^ — > 6*^, we have T/Ti = 

(1+0") ■ 

In other words, when a2 — > ai, W ^ (sec Eq. Q). 
Then Ti = T2 and T2 = Ti in this limit. In this case, 
the effective temperature is simplified to 



T 



+ 



(25) 



Finally, we make the following observations: 

i) it is interesting to note that the effective temperature 
is a weighted average of the subsystem temperatures after 
work extraction. 
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ii) the overall temperature is the same as the subsystem 
temperature when the latter are also equal to each other. 
This corresponds to v — Vd, which implies the well- 
known Curzon-Ahlborn efficiency 22 1. 

iii) At the global maximum of work, the conditions 
§^ = and = determine optimal values at and 

oai oa2 ^ 1 z 

as well as the condition 

^ = (26) 
aai aa2 



holds. So we have i^* = y^0/(,*, where now ^* is deter- 
mined from using optimal values and . Then it 
follows from Eq. ^ that 



2Ti T2 
Ti + T2 
This may be expressed as 

T*^ 2VF 



(27) 



(28) 



(1 + e 

iv) Fig. 1 shows the behaviour of subsystem and effec- 
tive temperatures as function of efhciency. Particularly, 
T shows a nonmonotonic trend. For a given value of rj, T 
in Fig. 1 corresponds to that engine set-up which yields 
the maximum work. From numerical calculations, it is 
observed that the temperature has a minimum at an effi- 
ciency which is bounded from below by Curzon-Ahlborn 
value. 
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FIG. 1: For T2 = 1 and Ti = 9, the subsystem temperatures 
Ti , T2 and the effective temperature T of the composite sys- 
tem evaluated at maximum work corresponding to a given 
efficiency. T is given by a weighted average over the subsys- 
tem temperatures and so its curve lies in between the curves 
for subsystem temperatures. All the three temperatures are 
equal at Curzon-Ahlborn efficiency. 



ferent definitions. We compare with our defintion two 
other definitions of the effective temperature from liter- 
ature, that are relevant to our system. The first can- 
didate is the spectral temperature [T^. This definition 
depends only on the energy probability distribution and 
the energy spectrum of the system and is applicable even 
for nonequilibrium situations. Thus for a non-degenerate 
spectrum, the inverse of spectral temperature is defined 
to be 



I 



= - 1- 



M 



-1 M 



E 



P,+P,^i \ l nP,-lnP,_i 

Ei — Ei^^ 



(29) 

where Boltzmann's constant has been set to unity. 
is the probability to occupy a level with energy Ei and 
index for the levels ranges from (ground state) to M . 
For our case of two TLS in the final state after work 
extraction, using the values {Ei\ = {0, a2,ai,ai -f 02} 
and {Pi} = {risi, r2Si, riS2, ?'2S2}, the (inverse) spectral 
temperature is explicitly given by 



1 



1 (i^- 



1 



I 



T[{v~v'^){l + x) T'^{\ + xy 



(30) 



where x = (ri + si — 2risi). The special cases include: 
Carnot limit, when v —t 6 and so P,. = (1 + x)T2] when 
V ^ Q,Ts ^ Q. Finally, for CA efficiency {u = ^/O), the 
spectral temperature is equal to the subsystems' temper- 
ature. In general, the behaviour of P, as shown in Fig. 2 
is quite different from the proposed definition. 




0.6 0.7 0.8 0.9 



FIG. 2: Comparison between other definitions of nonequilib- 
rium temperatures for r2 = 1 and Ti = 9. T is same as in 
Fig.l while Ts denotes the spectral temperature and Tc, the 
contact temperature, as defined in section IV. Inset shows the 
enlarged region around the point where all three temperatures 
are equal, which is at = 1 — \/9 = 0.6667. 



IV. EFFECTIVE TEMPERATURES: A 
COMPARISON 

The definition of nonequilibrium temperature is not 
unique for a given situation and one can envisage dif- 



The second definition we consider is also called as the 
contact temperature (Pc). If a general nonequilibrium 
system whose different parts may be at different local 
temperatures, is brought in contact with such a heat 
bath, that some parts of the system give heat to the 
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latter and some absorb heat from it, so that the net heat 
transferred between the system and the bath is zero, then 
the temperature of that bath defines [IS]- In other 
words, energy conservation holds for the system and dif- 
ferent parts of it come to a common temperature equal 
to that of the heat bath. Thus for two TLS, we impose 
that the total mean energy calculated with canonical dis- 
tributions for each TLS, corresponding to a temperature 
Tc, is equal to the final mean energy U = aiS2 + a2r2- 
The temperature obtained numerically is depicted in Fig. 
2. The behaviour of the contact temperature is closer to 
the proposed definition in regions where subsystem tem- 
peratures are equal or nearly to each other. However, 
towards the extreme values of the engine efhciency, the 
two temperatures take on different values. 

V. CONCLUSIONS 

The notion of temperature is well understood in the 
domain of equilibrium thermodynamics. However, its ex- 
tension to noncquilibrium situations is non-trivial. See 
for example [23 | for a review of effective temperatures 
in noncquilibrium situations. In this paper, we have dis- 
cussed a quantum heat engine in which two TLS prepared 
in different thermal states, undergo a unitary thermally 
isolated process and deliver work to an external work 
source. The final state of the two- TLS system is pas- 
sive (i.e. no further work can be extracted from it) but 



a noncquilibrium state where each subsystem may have 
a different local temperature. We have proposed a ther- 
modynamic definition to calculate effective temperature 
of the composite system in its final state. The obtained 
formula is very similar to the one expected on thermo- 
dynamic grounds. The proposed definition is compared 
with the spectral temperature, which seems to have a 
widely different behaviour. The other definition called 
contact temperature appears to have some semblence to 
our definition. All the three definitions converge for mu- 
tual equilibrium, but at Carnot limit or the vanishing ef- 
ficiency they differ from each other significantly. Future 
experiments on measurement of temperatures in such 
systems may decide between the different definitions. Fi- 
nally, it will be interesting to extend these ideas to more 
elaborate models such as involving entanglement between 
the TLS [2^. It is hoped that the present analysis will 
help to understand thermodynamic behaviour revealed 
by quantum heat engines. 
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